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By taking account of the alternation of structural parameters, we study bunching of impermeable 
steps induced by drift of adatoms on a vicinal face of Si(OOl). With the alternation of diffusion 
coefficient, the step bunching occurs irrespective of the direction of the drift if the step distance is 
large. Like the bunching of permeable steps, the type of large terraces is determined by the drift 
direction. With step-down drift, step bunches grows faster than those with step-up drift. The ratio 
of the growth rates is larger than the ratio of the diffusion coefficients. Evaporation of adatoms, 
which does not cause the step bunching, decreases the difference. If only the alternation of kinetic 
coefficient is taken into account, the step bunching occurs with step-down drift. In an early stage, 
the initial fluctuation of the step distance determines the type of large terraces, but in a late stage, 
the type of large terraces is opposite to the case of alternating diffusion coefficient. 

PACS numbers: 81.10.Aj, 05.70.Ln, 47.20.Hw, 68.35.Fx 



I. INTRODUCTION 

When a Si(OOl) vicinal face is tilted in the (110) direc- 
tion, two types of terraces, 1x2 and 2x1 terraces, appear 
alternately. In the 1x2 terrace, which we call Ta, the 
surface diffusion perpendicular to the steps is faster than 
that parallel to the steps. In the 2x1 terrace, which we 
call Tb, the relation is opposite. 

These terraces are separated by single steps. The type 
of the steps also changes alternately: the step at the 
upper side edge of Tb , which we call Sb , is rougher than 
that of Ta, which we call Sa- Then, parameters like 
kinetic coefficient and the step stiffness may change for 
the two kinds of steps. 

When a Si(001) vicinal face is heated by direct electric 
current, the vicinal face is unstable and step bunching 
occurs irrespective of the current direction P, Q. The 
type of large terraces between bunches is Ta with step- 
down current and Tb with step-up current. Cause of the 
step bunching is considered to be the drift of adatoms Q 
induced by the current. By using a one-dimensional step 
model, where the alternation of diffusion coefficient and 
that of kinetic coefficient are taken into account, Stoy- 
anov theoretically studied the stability of a vicinal 
face for pairing of steps. With large kinetic coefficients, 
the step pairing occurs irrespective of the direction of the 
drift, and the type of large terraces between step pairs is 
determined by the current direction. To study behaviors 
of step pairs, Natori and co-workers carried out numerical 
simulation of a similar one-dimensional step model 0, . 
With step-down drift, the step bunching occurs via co- 
alescence of step pairs, but the step bunching does not 



occur with step-up drift, which disagrees with the exper- 
iments [JQ. 

Recently, we carried out Monte Carlo simulation 0], 
where the alternation of anisotropic surface diffusion was 
taken into account. In contrast to the previous stud- 
ies 0, H|, tne s t e P bunching occurs irrespective of the 
direction of the drift as in the experiment 0- In 
our model , the steps are perfectly permeable, and the 
alternation of kinetic coefficient and the evaporation of 
adatoms are neglected. On the other hand, the steps are 
impermeable and both the alternation of kinetic coeffi- 
cient and the evaporation are taken into account in the 
previous models 0, El ■ Since there are many differences 
between the models, it is not clear what is the most im- 
portant factor to cause the different results. 

In this paper, we use a one-dimensional model of im- 
permeable steps to study the drift-induced step bunch- 
ing. In Sec. [HJ we introduce the model. In Sec. IIIII we 
analyze the model with alternating diffusion coefficient 
and perform a numerical simulation. We compare the 
results with our previous study 0. In Sec. lIVI we study 
the model with the alternating kinetic coefficient. The 
effect of evaporation, which is neglected in Ref. 0, is also 
studied in Sees. Illll and IIVI In Sec.^we summarize the 
results and give a brief discussion. 



II. MODEL 

We use a one-dimensional step flow model [J 0-0 . The 
y-coordinate is taken in the step-down direction. When 
the drift of adatoms is parallel to the y-axis, the diffusion 
equation of adatom density c(y, i) is given by 
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where D m is the diffusion coefficient in the mth terrace, F 
the force to cause the drift and r the lifetime of adatoms 
for evaporation. 

Boundary conditions at the mth step are given by 0,0 
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where c m is the equilibrium adatom density, K m ki- 
netic coefficient of the mth step and +(— ) indicates the 
lower (upper) side of the step. Ta and Tb appear al- 
ternately, and we assume steps with odd numbers are 
5b steps and those with even numbers are 5a steps 
(Figure The diffusion coefficients and the kinetic 
coefficients are given by (Di n ,Km) = (F>a,Ka) and 
(D 2n-i, K 2n -i) = (F>b,K b ). Bearing a Si(OOf) vicinal 
face in mind, we set Da > Db and Kb > Ka- The 
second terms in the right hand side of eqs. J2J) and © 
represent the adatom current through the step without 
solidification. With P m — > oo, the difference between 
c|, and c|_ vanishes and the step is called perfectly per- 
meable. With P m = 0, the step is called impermeable 
and the surface diffusion field is separated by the step. 
Hereafter we deal with the impermeable steps. 
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FIG. 1: A restructed Si(001) vicinal face. Short lines repre- 
sent dimers. 

When the neighboring steps interact with the potential 
£ TO ) the equilibrium adatom density at the mth step c m 
is given by 0, [Hj 
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where c° q is the equilibrium adatom density of an isolated 
step, fl the atomic area, y rn the position of the mth step 
and l m = (jjm+i — Dm) is the width of the mth terrace. 
If the step interaction potential C m is given by Q m = 



— A(\nl m _i + ln/ m ) as in a Si(OOI) vicinal face the 
exponent is v = f and A„ = SlA/k^T. 

By solving eq. with the boundary conditions, 
eqs. (J2J and @, in a quasi-static approximation (dc/dt = 
0), we determine the adatom density. The velocity of the 
mth step is given by 
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In the following, to see how the step bunching changes 
with the alternation of parameters, we separately study 
the effect of diffusion coefficients and that with kinetic 
coefficients. 



III. STEP BUNCHING WITH ALTERNATION 
OF THE DIFFUSION COEFFICIENT 

We use the model of impermeable steps to study the 
step bunching with alternation of diffusion coefficient. 
We compare results wtih the previous study to see the 
effect of step permeability. For simplicity, the alternation 
of kinetic coefficient is neglected. 



A. Step bunching induced by the drift in a 
conserved system 

We first study the step bunching in a conserved system. 
Without the evaporation, the step velocity is given by 
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where / = F/k B T. 

In a vicinal face with the step distance I, the step ve- 
locities are given by 
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Since the diffusion coefficient Da is larger than Db, 5a 
steps advance and 5b steps recede with step-down drift 
(/ > 0). With step- up drift (/ < 0), the direction of the 
step motion is reversed. 
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In the initial stage of the instability, formation of large 
Ta(Tb) terrace with step-down (step- up) drift is ex- 
pected due to the step pairing. To study the motion 
of pairs of steps, we carry out numerical integration of 
cq. (JBJ. In addition to the step distance I, there are 
two characteristic lengths in the vicinal face. One is 
f^ 1 ~ k^T/eE [l^, which is a characteristic length de- 
termined by the drift and usually much larger than the 
step distance. The other is Da/K, which is character- 
istic length determined by step kinetics. When the step 
distance I is much smaller than Da/K, the step kinetics 
is more important than the surface diffusion. We assume 
that / _1 is larger than Da/K and I, and study the time 
evolution in two cases; Da/K < I and I < Da/K. 

300 n — ■ 1 t i ■ 1 




FIG. 2: Time evolution of step position. The step distance 
I is larger than Da/K and Db/K. The drift direction is (a) 
step-down drift and (b) step-up, and the drift velocities satisfy 
\fl\ = 0.1. The number of steps is 32 and the system size is 
16 with the periodic boundary condition. An = 5 x 1CF 3 with 
v = l, Da/KI = 0.2, D B /Kl = 0.02 and I = Qc° q t. 

Figure[5]represents the time evolution of step positions 
with Da/K <C I <C This is the case in which the 

diffusion is slow and controls the time evolution. Solid 
lines represent the evolution of the positions of Sb and 



dotted lines represent Sa- The characteristic length / _1 
is / _1 = 200 and the scaled time t is given by t = OCg q i. 
The initial step distance is about the same, but with a 
random fluctuation. In the initial stage, the step pairing 
occurs irrespective of the drift direction. Large terraces 
are Ta with step-down drift (Fig. [2Ja)) an d Tb with 
step-up drift (Fig. [^bThwhich agrees with eq. iJJJ and 
the previous analysis 0,1!, El- 

In a late stage, large bunches appear irrespective of the 
drift direction. The type of large terraces is the same as in 
the initial stage. Free single steps or pairs do not exist on 
the large terraces and the bunches grow via coalescence 
of small bunches. Though the formation process of large 
bunches is similar in both cases, the step density with 
step-down drift is higher than that with step- up drift, and 
the growth is much faster. The ratio of the growth rate 
of step bunches is larger than the ratio of the diffusion 
coefficients. In the numerical study of Natori et al. 0,0, 
the formation of large bunches with step-up drift did not 
occur in contrast to our simulation. For the very slow 
growth rate of step bunches, their simulation time might 
be too short to produce large bunches. 

Figure|3]represents the step bunching with I <C Db/K, 
i.e., the step kinctics-contral case. The amplitude of the 
initial fluctuation is the same as that in Fig. [5] When 
the drift is in the step-down direction (Fig. EJa)), the 
equidistant step train is unstable and the step bunch- 
ing occurs. The terrace type between step bunches is 
determined by the initial fluctuation and both types of 
large terraces coexist. When the drift is in the step-up 
direction (Fig.JSfb)), the step bunching does not seem to 
occur. The results are very different from Fig. [5] 

To find the reason that the step behavior changes with 
the kinetic coefficients, we analyse the change of alter- 
nating terrace width. We assume the width of Ta is 
l 2n = 1 + Al/2 and that of T B is hn-i = I - Al/2. Fig- 
ure0|rcprcscnts the time derivative of Al for large kinetic 
coefficients. Parameters are the same as in Fig. [21 With 
step-down drift (Fig.^Ja)), the initial vicinal face Al = 
is unstable and Al increases up to the stable fixed point 
near Al = 21. The surface consists of large Ta terraces 
and small Tb terraces, i.e., tight step pairs. With step- 
up drift (Fig. Bib)), the initial stage is also unstable and 
Al decreases to the stable fixed point near Al = —21. 
The surface consists of small Ta terraces and large Tb 
terraces. The results agree with the initial pairing stage 
of Fig. H 

Figure [3] represents the time derivative of Al for small 
kinetic coefficients. With step-down drift (Fig. [5ja)), 
there are three fixed points. The fixed point (open cir- 
cle) near the center is unstable and two other fixed points 
(filled circles) are stable. Since the unstable fixed point 
is very close to Al = 0, the final stage is not unique if 
the initial fluctuation is included. In Fig.[3{a), the initial 
random fluctuation is not negligible and some terraces 
move to the fixed point with positive Al and two types 
of terraces coexist. 

With step- up drift, there is only one stable fixed point 
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FIG. 3: Time evolution of step position with step distance 
I <C Da/K, Db/K with (a) step-down drift and (b) step- up 
drift. The number of steps is 32 and the system size is 16 
with the periodic boundary condition. A,y = 4 x 10~ 3 with 
v = l, |//|=0.1, Da/KI = 200 and D B /Kl = 20. 




FIG. 4: Time derivative of the terrace width change Al with 
(a) step-down drift and (b) step-up drift. Parameters are the 
same as in Fig. [2] 



with a small negative Al. From the condition V m in 
cq. ©, the difference of step distance Al in the fixed 
point near Al = is given by 



M = Kl(D B - D A ) 
I 2D A D B ' 



(9) 



where we neglected the step-step repulsive interaction. 
The sign of Al is determined by the diffusion coefficients 
and independent of the drift direction as seen in Fig.^b). 



B. Step bunching induced by the evaporation 

Since the experiments Q, Q were carried out at high 
temperatures, the evaporation of adatoms may not be 
negligible for the step bunching. 

We first neglect the drift of adatoms and see if the step 
bunching occurs with the evaporation. For simplicity, we 
consider the limit of the fast step kinetics, K — > oo. If 



the step repulsion is absent, the step velocity is given by 
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where x m = V ' D m r is the surface diffusion length in the 
mth terrace. With the same step distance I, the step 
velocity is given by 
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where xa(^b) represents the surface diffusion length in 
Ta(Tb) and the equidistant train of steps is a steady 
state. When the terrace widths change alternately and 
are given by hn = I + Al/2 and fen-i = I — Al/2, from 
cq. IjlOU the time evolution of Al is given by 
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FIG. 5: Time derivative of the terrace width change Al with 
(a) step-down drift and (b) step-up drift. Parameters are the 
same as in Fig.|H] 



The vicinal face is marginal to the perturbation. If the 
repulsive interaction is taken into account, the vicinal 
face is stable and the step pairing does not occur. With 
the alternation of diffusion coefficient, the evaporation 
alone does not cause the step bunching. 

When both the drift and the evaporation are present, 
the step velocity is given by 

Vm (An-l - Dm)fc m 

n 2 



D m -ia m -i(c m cosha m _iZ m _i - e 



-/Jm-l/2 



Cm-l) 



sinha TO _iZ m _i 

-DfnCXm (c m COSll CX 7n l m 6 ^ m Cm+l)? 



The parameter a m is defined by 



(13) 



(14) 



o^ 1 is the characteristic length of the diffusion field in 
the mth terrace. 
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FIG. 6: Time evolution of step position with evaporation 
with (a) step-down drift and (b) step-up drift. The number 
of steps is 16 and system size is 8. The diffusion coefficients 
are Db = 2 and Da = 10. Other parameters are / = 0.1 
v€lA/knT = 10" 4 with v = 1, and r = 50. 



We carry out numerical integration of cq. I|13f) . The 
time evolution of step positions (Figure |SJ) shows the 
drift-induced step bunching with the evaporation, where 
the characteristic length scale in Ta is a^ 1 = 14.9 and 



that in Tb is a 2 



.9. Irrespective of the drift di- 



rection, very fast pairing of receding steps occurs at the 
initial stage. The type of large terraces is the same as in 
Fig and not changed by the evaporation. 

Large bunches appear by coalescence of small bunches. 
With step-down drift, coalescence of step bunches oc- 
curs successively, but the time interval of the coales- 
cence increases with increasing the terrace width between 
bunches (Fig. Efa)). With step-up drift, coalescence of 
step pairs does not occur until t « 150. However, once 
the coalescence of step pairs starts, the interval does not 
seem to increase much (Fig. Efb)). The difference of the 
growth rate of bunches is much smaller than that without 
the evaporation. 
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IV. STEP BUNCHING WITH ALTERNATION 
OF KINETIC COEFFICIENT 

In a Si(OOl) vicinal face, not only diffusion coefficients 
but also the type of steps changes alternately. In this sec- 
tion, we study the step bunching due to the alternation 
of kinetic coefficient and compare the result with Scc. lIIII 



A. Step bunching induced by the drift 

To focus on the alternation of step kinetics, we neglect 
the alternation of diffusion coefficients. First, we neglect 
the evaporation. Without the evaporation, the step ve- 
locity is given by 



100 
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where D s is the diffusion coefficient. 

In a vicinal face with equidistant steps, the step veloc- 
ities are given by 



V 2n = -V 2n - 



nD s K A K B c°f 2 (ef l 



\f(K B -K A ) 



,(16) 

ffAS'B 

where c/a and <?b are 

g A = [K A (DJ -K B )+K B (D s f + K A )e n ], 

g B = [K B (D s f-K A )+K A (D s f + K B )e fl ]. (17) 

Since the kinetic coefficient K B is larger than K A , Sa 
steps advance and Sb steps recede. Alternation of large 
Ta terraces and small Tb terraces is expected by forma- 
tion of step pairs. 

When the kinetic coefficient K B is large, D s /Kb <C I, 
time evolution of step positions is shown in Fig. [7J With 
step-down drift (Fig. Ufa)), most steps form pairs and 
large Ta terraces appear in the initial stage, but forma- 
tion of triplets also occurs and a large Tb is produced. 
In a late stage, all large terraces become Ta after coales- 
cence of small bunches. With step-up drift (Fig. Efb)), 
the pairing does not seem to occur. 

Figure |H| represents the growth rate of Al, which is 
defined in the same way as in Sec. lIII Al With step-down 
drift, an unstable fixed point with a small positive Al and 
two stable fixed points with a large amplitude of |AZ| are 
present. From the condition V m = in cq. (|I5f) . the 
difference of terrace width at the fixed point near Al = 
is given by 



B 



Al _K A - K 
T ~ K A + K B 



J/- 



(18) 



where we have used /!< 1 and neglected the step-step 
repulsion. In cq. 10 Al is determined only by structural 
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FIG. 7: Time evolution of step position with alternating 
change of kinetic coefficient with step distance D s /Ka <C 1. 
Drift is in (a) step-down direction and (b) the step-up di- 
rection. The number of steps is 32 and the system size is 
32 with the periodic boundary condition. Parameters are 
M = 2 x 1(T 2 with 5 = 1, \fl\ = 0.2, D s /KaI = 1.0 x 10~ 2 
and D s /KbI = 1.0 x 10~ 3 . 



parameters and independent of the drift, but in eq. 1181) 
Al depends on the drift and changes the sign with the 
drift direction. Expect for the sign of Al in the unstable 
fixed points, the form of dAl/dt is the same as in Fig. [S] 
The initial vicinal face moves to the steady state with the 
large negative Al if the initial fluctuation is small, but it 
can move to the other steady state with large Al if the 
initial fluctuation is large. In Fig. [3a), formation of a 
large Tb in the initial stage is due to a large fluctuation 
of terrace width. 

With step-up drift, there is only one fixed point with 
a negative Al. The steady state is stable. Since the 
difference of terrace width Al is small, the surface looks 
like the original vicinal face in Fig.0h). 

With the alternation of diffusion coefficient, the form 
of dAl jdt changes with decreasing the kinetic coefficients 
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FIG. 8: Time derivative of the terrace width change with 
(a) step-down drift and (b) step-up drift. Parameters are the 
same as in Fig. [7| 



(Figs. |H and On the other hand, with the alternation 
of kinetic coefficient, the form of dAl/dt does not change 
even if D s /K changes. Thus the bunching behavior is 
insensitive to the ratio D s /Kl. 



B. Step bunching induced by the evaporation 

If the evaporation is present and the drift is absent, 
the step velocity is given by 



(A m +i svahl m /x s + coshJ m /x s ) 



'n+1 



h m +l(lm) 
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where A TO = D s /K rn x s , and h m (l) is 

h m (l) = (A m _iA m + l)sinh//x s 

+ (A m + A r „_i)coshZ/a; s . (20) 



Due to the evaporation, both Sa steps and Sb steps re- 
cede. When the steps are equidistant, the difference of 
step velocities is given by 



V 2 n - V, 



2n-l 



( Aa — Ab ) tanh I / x s 



(1 + A A A B ) tanhZ/a; s + (A A + A B ) 



(21) 



Since Aa is smaller than A B , 5b recedes faster than Sa- 
The vicinal face is unstable, and large Tb terraces and 
small Ta terraces appear by step pairing. 
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FIG. 9: Time evolution of step position. a- s = 15, Ka = 
0.125, K B = 1.0, D B = 1.0. The number of step is 16 and 
the system size is 8 with the periodic boundary condition. 
Ac = xl0~ 4 with v = l. 

Figure[5]represents the time evolution of step positions, 
which is obtained by numerical integration of cq. i|19fl . 
In the simulation, the difference of kinetic coefficients is 
very large. In very early stage, pairing of steps occurs 
and large Tb terraces appear. Via coalescence of step 
pairs, step bunches appear. The type of large terrace is 
the same as that in the step bunching induced by the 
drift. In large terraces, isolated step pairs are present. 
When a collision between a step bunch and a step pair 
occurs, another step pair leaves from the upper side of 
the step bunch. A step pair breaks into single steps, but 
only temporarily. Repea ting of such collisions is seen in 
other bucnhing systems |l4lll5l|. Here a step pair is the 
fundamental unit. 



V. SUMMARY AND DISCUSSION 

In this paper, we studied the drift-induced step bunch- 
ing with the alternation of structural parameters: the 
diffusion coefficient and the kinetic coefficient. 

With the alternation of diffusion coefficients on con- 
secutive terraces, the step bunching occurs irrespective 
of the drift direction if the kinetic coefficients are large. 
The type of large terraces is determined by the drift di- 
rection. When the kinetic coefficients are small, the step 
bunching occurs with step-down drift and the initial fluc- 
tuation of step distance influences the type of step pairs. 
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The growth rate of bunches is much faster with the step- 
down drift, but the difference of growth rates decreases 
with the evaporation. Without drift of adatoms, the 
evaporation does not induce bunching. 

With the alternation of kinetic coefficient at consecu- 
tive steps, the step bunching occurs with the step-down 
drift and does not occur with the step-up drift. The type 
of large terraces is influenced by the initial fluctuation 
of step distance. The evaporation induces step bunching 
even if the drift is absent, in contrast to the diffusion 
coefficient case. 

In the experiments P, Q , the type of terraces between 
bunches changes with the drift direction when the initial 
step distance is large. From our result we may conclude 
that the alternation of diffusion coefficient is essential 
to the bunching and that of kinetic coefficient is not so 
important. Also the kinetic coefficients are large if the 
steps are impermeable. Impermeable steps with large 
kinetic coefficients are effectively equivalent to the per- 
meable steps However, when the step distance is 
smaller than the critical value, the step bunching occurs 
only with the step-up current in the experiment [l( , which 



is not explained by our model. In the previous stud- 
ies [13, Ell , the dift direction to cause the step bunching 
on a Si(lll) vicinal face is affected by the step perme- 
ability. In a similar way, the disagreement in a Si(OOl) 
vicinal face may be explained by the step permeability, 
which remains to be explored. 

The difference of growth rate of step bunches as the 
change of current direction is very small in the experi- 
ment 0]. The evaporation may play an important role, 
but the difference does not vanish only with the evapo- 
ration. In the Monte Carlo simulation pfl, in which steps 
are permeable, a similar behavior has been observed. We 
still do not undestand the origin of such a large difference. 
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